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Abstract 

The number of absorption and Raman optical modes for each Wyckoff orbit in the high- 
temperature tetragonal (I4/mmm) parent lattice of the Ruddlesden-Popper compounds CasX207, 
with X=Mn or X=Ti is given. We analyze the effect of sequential perturbations which lower the 
symmetry to Cmcm and Cmc2i and finally include magnetic ordering. We determine the power 
law behavior (within mean field theory) for the cross section for photon absorption and Raman 
scattering of modes which appear as the symmetry is successively lowered. In the Cmc2i phase we 
give a symmetry analysis to discuss the magnon-phonon coupling which in other systems gives rise 
to "electromagnons." From our results we suggest several experiments to clarify the phase diagram 
and other properties of these systems. 

PACS numbers: 61.50.Ks,61.66.-f,63.20.-e,76.50.+g 
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I. INTRODUCTION 



The Ruddlesden-Popper (RP) compounds [1] are those of the type A n+1 B n C 3n+ i, where 
the valences of the ions are A = +2, B = +4, and C = —2 (usually C is oxygen). At 
high temperature the crystal structure for many of these compounds is tetragonal, but the 
structure can distort in several ways as the temperature passes through structural phase 
transitions. In the lowest symmetry phase, these systems are usually ferroelectric and can 
also be magnetically ordered, although, in contrast to many magnetoelectrics, {2-7] these two 
degrees of freedom do not appear at the same temperature. Nevertheless the RP systems 
do have many interesting couplings between their magnetic and dielectric properties. For 
example, their spontaneous polarization can be rotated by the application of a magnetic 
field and the weak ferromagnetic moment can be rotated by the application of an electric 
field, js, 9] The symmetry of the RP's is very similar to that of the Aurivillius compound 
SrBi2Ta20g analyzed in detail by Perez-Mato et al. 10[ and many of our results should also 
apply to it. Although there have been several first principles calculations of the properties of 
CMO and/or CTO, the calculations of CF are most relevant to the present paper. 



FIG. 1: (Color online) The tetragonal phase of CMO or CTO (space group I4/mmm). The oxygen 
ions are the blue spheres, the Ca ions are the open green circles, and the Mn or Ti ions (not shown) 
are at the centers of the octahedra of oxygen ions. There are two molecules per conventional unit 
cell. 
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Here we will be interested in the two RP compounds Ca 3 Mn 2 7 (CMO) and Ca 3 Ti 2 07 
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(CTO), which consist of bilayers of oxygen octahedra enclo sing a M n or Ti ion. It is believed 
that the high-temperature phase is tetragonal (I4/mmm) 13lJl4j and the low-temperature 
phase is orthorhombic (Cmc2i) 15[ and in Fig. [2] we show the three possible paths that con- 
nect these structures. [8] First principles calculations jsj indicate that the most likely scenario 
involves the appearance, as the temperature is lowered, of the intermediate state Cmcm at a 
temperature we denote T>. In fact, the Cmcm phase has been observed in the isostructural 
compounds LaCa 2 Mn 2 07[ii]] at and below room temperature and Bio.44Ca 2 . 56 Mn 2 07[l7| at 
room teperature. Accordingly, we will assume that this scenario applies to CMO and CTO 
although for these systems no intermediate phase has been observed up to now. However, re- 
cent measurements on ceramic CMO find a clear pyroelectric signal consistent with the onset 
of ferroelectric order close to T> = 280K. 18] This transition is identified as the temperature 
at which Cmc2! appears. Since this ferroelectric transition seems to be a continuous and 
well-developed one and since a direct continuous transition between I4/mmm and Cmc2! is 
inconsistent with Landau theory, 19] the seemingly inescapable conclusion is that the phase 
for T slightly greater than T> is not I4/mmm, but is some phase which does not allow a 
spontaneous polarization. Thus the phase at temperature just above T > may be the Cmcm 
phase. In this scenario the Cmc2i phase would appear at a lower temperature T< (which 
experimental data implies is quite close to T>). The purpose of the present paper is to 
discuss the symmetry of the zone center phonons and the Raman scattering cross section 
as they are modified by the successive lowering of symmetry as the temperature is reduced 
through the structural phase transitions at T> and T<. Such a symmetry analysis can be 
done using standard group theory methods. Here we will supplement that analysis with an 
analysis of the phonon cross sections from which we determine their power law dependence 
on the order parameters (OP's) which characterize the distortions from tetragonal symmetry 
to the lower symmetry structures. This enables us to predict the temperature dependence 
of the newly allowed cross sections which appear at the structural phase transitions. We 
then give a similar analysis of the magnon-phonon interaction. As a result of this coupling, 
the magnon absorption cross section (which usually results from a magnetic dipole allowed 
transition) is enhanced by now being electric-dipole allowed. This enhancement has led some 
authors 20M25| to call such magnons "electromagnons." 

Briefly, this paper is organized as follows. In Sec. II we give a symmetry analysis 
of the zone center phonons in the tetragonal (I4/mmm) phase and in Sees. Ill and IV 
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we give similar analyses for the orthorhombic Cmcm and Cmc2 1 phases. We discuss the 
photon absorption and Raman scattering of these modes. For modes which only appear as 
a consequence of a phase transition we also determine the power-law dependence of their 
optical cross sections on the emergent OP's. In Sec. V we discuss the coupling of these 
modes to magnetic excitations (magnons). Our results are summarized in Sec. VI. 

I4/mmm(139) 



X, 



Cmca (64) 



F2mm (42) 



X 



Cmcm (63) 



A2 1 am (36) 



FIG. 2: The group-subgroup structure of CMO and CTO. We indicate the irreducible representa- 
tions (irreps) by which the tetragonal phase is distorted to reach one of the intermediate phases. 



The numbering of the space groups follows Ref. 



26l . The superscripts on the irreps give the parity 



with respect to spatial inversion. The irrep is at zero wave vector and the X irreps are at 
the wave vectors q = (1/2, ±1/2, 0). In the final step to reach the A2iam (equivalent to Cmc2i) 
structure the other two irreps are introduced simultaneously, |9] so that Cmc2i can be regarded as 
the result of simultaneous condensation of the three irreps. 



II. SYMMETRY ANALYSIS FOR THE TETRAGONAL STRUCTURE 

The occupied sites in space group I4/mmm for the n = 2 RP system are given in the 
first three columns of Table IIII[ below. The character table for the little group of q = for 
the I4/mmm tetragonal structure is given in Table [H For the one dimensional irreducible 
representations (irreps) the characters are the matrices. For the two dimensional irreps and 
for the irreps for the X points at q = (1/2, ±1/2, 0) (in rlu's) the matrices are given in Table 



ED 28] 
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TABLE I: Character table for the little group of q = for the I4/mmm tetragonal structure. The 



operator numbers are those of Ref. l27|. The bottom row gives the dipole moment vector P or the 
Raman tensor R, where a, b, c, etc. are arbitrary constants. For the two dimensional irreps the 
two entries are the values of P and R for each of the two degenerate states. The entry is left blank 
for irreps which are optically inactive. 
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We now discuss briefly the introduction of order parameters (OP's). Starting from the 
high-temperature tetragonal, I4/mmm, structure, distortions of the irreps , X% and X% 
are condensed as the temperature is lowered to eventually reach the low-symmetry low- 
temperature Cmc2! structure. A tabulation of these symmetry-adapted distortion vectors, 
, for irrep X is given in Table IHIi In the dominant component of the distortion of irrep 
X£ the octahedra are rotated about the crystal c axis like interlocking gears. The rotations 
in the upper (z = r) and lower (z = — r) planes of a bilayer are equal, so in the notation 



of Refs. 



29 and 



301 these rotations are both denoted (0,0,©). The dominant component 
of the distortion of irrep X% consists of 1) a rotation about the tetragonal [1,0,0] with the 
rotation alternating in sign as one moves along the [1,0,0] direction and 2) a rotation about 
the tetragonal [0,1,0] direction with the rotation alternating in sign as one moves along the 
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TABLE II: Representation matrices M^ G \0) for the operators O which are the generators of the 
irrep G. Here r' = Or and the cr's are the Pauli matrices. For the irreps G at zero wave vector the 
matrices for the translations Ti = (x + 1, y, z), T2 = (x, y + 1, z), and T3 = (x + ^, y + ^, z — ^) 
are unity. For the irreps at the X wave vectors the matrices are M(Ti) = M(T2) = — 1 and 
M(T3) = — o z . The first row and column for the X irreps refer to wave vector qi These matrices 
are related to those of Ref. |l9| by a unitary transformation which, for zero wave vector takes 
(<t x , cr y , cr z ) into (— cr x , —cr y , a z ) and for the wave vector X takes (cr x , a y , a z ) into {—a z ,a y , —cr x ). 
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[0,1,0] direction. Since this distortion is an even function of z, we identify this distortion as 
($, $, 0) in the upper layer and (— $, — $, 0) in the lower layer. 

If we assume that the distortions are small, then the crystal structure in any of the phases 
can be approximated as a distortion from the tetragonal structure which is written as linear 
combination of the distortions of the three irreps listed in Table IHIt 

$ = E Q%>*SP , 

m,X 

where ^ffl is the 12 component vector of irrep X as in Table IHIl Since the notation to 
distinguish the various irreps which we will encounter is somehwat complicated, we summa- 
rize the notation in Table IIVI If we regard the basis vectors as fixed, then application 
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TABLE III: Basis functions and ^ for the distortion under irrep X, where X = 5,3,2 

indicates irrep Tg , X% , and X^", respectively. |28] For each site we give the three components of the 
vector displacement. We assume normalization so that the sum of the squares of the components 
is unity. The values of the structural parameters are given in Ref. 1 131 : p = 0.311, £ = 0.100, 
X = 0.205, and r = 0.087. Since the a and b sites have the same symmetry we will refer to them 
both as a sites. 
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of an operator O in the symmetry group of the tetragonal space group to $ will induce a 
transformation of the OP's as 

k 

where the representation matrices ~M.( X \0) for the all the relevant operators O can be 
constructed from those given for the generators of the irrep X given in Table [TTJ This 
equation defines the symmetry of the OP's. 



TABLE IV: Notation for order parameters. First column: full notation, subsequent column(s): 
abbreviated notation. 
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As explained in Ref. |9| these distortions form two families, one of which involves wave 
vector qx = (1/2, 1/2, 0) and the other of which involves wave vector q 2 = (1/2, —1/2, 0). 31 | 
The wave vectors (1/2, ±1/2,0) correspond to a doubling of the size of the primitive unit 
cell. The first family is described by the three OP's Q3 (qi), Q^(qi), and Z oc Q 1 5 +Q 2 5 , 
and the second by the OP's (q2), Qt (02); an d Y oc Q 1 5 — Q 2 5 • In the absence of any 
external fields, when the sample orders, one expects an equal population of domains of the 
two families. However, it is possible to select a single family by invoking the coupling jg] 



V = ae X j / [Q fc (qi) 2 - <5fc(q 2 ) 



(1) 



where e xy is a shear component of strain, and k = 2 or 3. One can utilize this interaction 
by cooling the sample through T> in the presences of a shear stress conjugate to e xy . This 
will preferentially select qi or q 2 , depending on the sign of the stress and the sign of the 
constant a. It is possible that the selection of the wave vector can also be accomplished by 
cooling in an electric or magnetic field. jo| 
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Our aim is describe the evolution of the cross section for a) zero wave vector photon 
absorption and b) Raman scattering in terms of the various OP's. 

The single photon absorption process results from the electromagnetic coupling Ve = 
— p • E, where E is the electric field of the incident electromagnetic wave and p is the dipole 
moment operator of the system. The symmetry indicated in Tables [J and [III indicates that 
in the tetragonal phase one has 

Px = a x Q 5 (l) , p y = a y Q 5 (2) , p z = a z Q^(0) , 

with a y = a x . Note: in reality we ought to index the modes by an additional index to 
distinguish between different modes of the same symmetry Thus, for instance, 

«(r 5 -) 

Px = a k,xQh{])k , (2) 
fc=i 

where n(T^) is the number of occurrences of , which, in this case, is [see Eq. (]4j), below] 
n(r^) = 7. For simplicity we will usually not display this index explicitly, but in all our 
results sums over such an index (or indices) are implied. We discuss the use of Eq. (J2J) 
in Appendix |A] The absorption cross section is proportional to the square of the scalar 
product of the polarization vector of the incident photon and the dipole moment vector for 



the irrep in question. [32|J33| (For the two dimensional irrep one has to sum the squares of 
the two scalar products.) 

Likewise, the Raman scattering process results from the electromagnetic coupling Ve = 
— H(3 1 0l MEpE 11 where a. is the polarizability tensor, which depends on the phonon dis- 
placements as 



(o) . \ - dap-y 



We consider the cross section (or intensity) for a process in which the incident photon has 
an electric field vector in the direction and the scattered photon has an electric field 
vector in the 7 direction (we call this the ^-7 polarization) and a phonon Q T is created 
(or destroyed). The cross section for this process is proportional to the square of the tensor 



product of t 



re Raman tensor da^/dQ T with the polarizations of the incoming and outgoing 



photons. [32lj33H (For the two dimensional irreps one sums these squares over the two states 
of the irreps.) For both absorption and scattering we implicitly assume that the photon 
wave length is very large in comparison to the size of the unit cell. 
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Now we determine which modes (at zero wave vector) are allowed in the tetragonal lattice. 
Using standard methods we find that the reducible representations T corresponding to the 
space of zero wave vector displacements of the Wyckoff orbits a, e, and g, are decomposed 
into their irreps as 

r(a) = r 5 - + r 3 - , 

r(e) = r 5 - + r3 +r+ + r+, 

v(g) = 217 + r 4 - + r- + 2r+ + r+ + r+ , (3) 

so that the reducible representation T corresponding to the space of the vector displacements 
of all 12 ions in the primitive unit cell has the decomposition into irreps as 

T = 2T(a) + 3T(e) + T{g) = 717 + T 4 " 

+er 3 +5r+ + 4r+ + r+ . (4) 

The allowed absorption modes [x and y from 17(0) and z from 17(0)] are listed in Table PVT. 
Thus we find that in the tetragonal phase there are 20 absorption modes, 6 from 17 which 
have z polarization and 7 pairs of doubly degenerate modes (one with x polarization and 
one with y polarization) coming from 17 . Note that of these 20 absorption modes three are 
acoustic (zero frequency) and are not observed in IR experiments. There are also 15 Raman 
active modes, five nondegenerate modes coming from irreps Tf and and 5 pairs of doubly 
degenerate modes coming from Tf. Because the lowest symmetry sites of the space group 
I4/mmm are not occupied in CMO or CTO, it happens that there are no modes for which 
the Raman polarization xy is nonzero. Indeed, this result can be used as a signature of the 
tetragonal phase of CTO or CMO. 

III. ORTHORHOMBIC SPACE GROUP CMCM (# 63) 
A. Analysis based on Cmcm symmetry 

We now consider the modifications in the optical spectrum when we condense irrep X% 
at wave vector qi = (1/2,1/2,0) to go into the space group Cmcm. The transformation 
from tetragonal to orthorhombic coordinates is illustrated in Fig. [3] and is given by 

X = zt, Y = — -, 
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TABLE V: Number of modes Q of the tetragonal phase with their polarizations x, y, or z, for 
Wyckoff orbits, RP(re) for n = a,b,c. The phonon modes Qf(l) and Qf(2) are degenerate in 
energy. To identify with modes in the orthorhombic phase zt — > Xo, (xt, yt) — > (Yo, Zo), where the 
subscript O (t) refers to the orthorhombic (tetragonal) coordinate system. Pol. denotes polarization 
in a notation where xy indicates that the Raman tensor is R = const(r/ + ji), XoYo indicates that 
R = const (IJ + J I), where lower case unit vectors are tetragonal and capitals are orthorhombic, 
and similarly for other polarizations. Also Yq — Zq indicates that R = const. ( J J — KK) and X , 
(Yq + Zq) indicates that R = all + b{JJ + KK), where a and b are constants. (The orthorhombic 
polarizations are those when the tetragonal symmetry is broken.) The row labeled TOT gives the 
total number of modes for 2 a orbits, 3 e orbits, and 1 g orbit. 
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a) Here and below the absorption modes include an acoustic mode for each polarization. 

xt + yt 1 

z o = 7 , (5) 



where here and below we use lower case letters for tetragonal coordinates and capitals for 
orthorhombic coordinates. The unit cell is twice as large as the t unit cell. That is why 
x t ± yt is divided by 2 rather than by a/2. In Table I VI I we give the character table for the 
little group of q = for Cmcm as well as the characters for the reducible representation 
associated with the vector displacement over each of the Wyckoff orbits. 

To find out how many Cmcm optically active modes there are, we need to convert the 
tetragonal coordinates into the orthorhombic coordinates but we do not need to explicitly 
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FIG. 3: (Color online) Tetragonal (large black arrows) and orthorhombic (small red arrows) axes. 

TABLE VI: Character table for the little group of q = for Cmcm and (in the last row) how they 
transform. The last three columns give the character vectors for displacements for the Wyckoff 
orbits of the RP lattice. The 1^ modes are Raman active with Raman tensor R = all+bJ J+cKK . 
All coordinates are orthorhombic. The modes of irrep Tg transform like XY Z and are silent. 
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include the distortion that makes the orthorhombic structure. From the parent tetragonal 
orbit a we have the orthorhombic orbit 

ri = (0,0,0) t = (0,-^ ~) , 
r 2 = (l,0,0) t =(0,^) o , 
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where the subscript t means tetragonal and the subscript o means Cmcm. We find that 

T(acmc m ) = r+ + r+ + r+ + 17 + 17 + 17 . (6) 

For the parent tetragonal orbit e we have 

ri = (0,0, z) t = (z,~,--) , 

. 1 1, 
r 2 = (1,0, z) t = {z,-,-)o , 

1 1. 

r 3 = (0,0,-z)t = (-z, --,--) , 

1 1. 

r 4 = (1,0, -z) t = {-z, -, -) , 

and we find that 

r(e Cmcm ) = 217 + r+ + 17 + 217 + 217 + 17 

+r 3 " + 2r 4 - . (7) 

Similarly, for the parent tetragonal orbit g we find that 

T(gcmcm) = 217 + 417 + 217 + 4r+ + 217 

+417 + 217 + 417 . (8) 

Thus the decomposition of the reducible representation for the displacements of all ions in 
the unit cell is 

man ) + 3r(e c in cm ) + ^{9 Cmcm) 
= 1017 + 717 + 717 + 1217 + 1017 
+917 + 517 + 1217 . 

So in the Cmcm phase we have 10 modes with Xq polarization, 12 modes with Yq polar- 
ization, and 9 modes with Zq polarization, as given in Table IVHI (The degeneracies of the 
tetragonal structure are removed in the Cmcm structure.) Results for the Raman-active 
modes are also given in Table I VIII 

B. Analysis based on lowering of tetragonal symmetry 



We now discuss the phonon modes by analyzing the perturbative effect of distortions from 
the tetragonal parent structure. Compared to the previous analysis, this analysis provides 
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TABLE VII: As Table [V] Modes Q in the orthorhombic Cmcm phase for the Wyckoff orbits 
descended from the tetragonal orbits. The irreps are those of Table IVII for the Cmcm phase. 
In the next to last row we give the total number of modes for 2 a orbits, 3 e orbits and 1 g orbit. 
In the last row we repeat the results for the total number of modes for the tetragonal system. 





Absorption 


Raman 


IRREP 


rr 


r 2 


r 4 " 


rf 


1 2 


1 3 


1 4 


Polarization 


Xo 


Zo 


Yo 


y2 v 2 V 2 


XoZo 






RP(a) 


1 


1 


1 


1 





1 


1 


RP(e) 


2 


1 


2 


2 


1 


1 


2 


RP(g) 


2 


4 


4 


2 


4 


2 


4 


TOT 


10 


9 


12 


10 


7 


7 


12 


TET 


6 


7 


7 


4 


5 


1 


5 



information on how the intensity of the new Cmcm modes depends on the amplitude of 
the orthorhombic distortion from tetragonal symmetry. Let us suppose that Q% (qi) has 
condensed due to the phase transition at temperature T = T> into the Cmcm phase. (Near 
this phase transition one can choose modes so that only a single mode of this symmetry 
is involved. As the temperature is lowered below T> other modes of this symmetry will 
progressively be introduced.) To see how the existence of the new order parameter can 
modify modes of the tetragonal system consider the Hamiltonian for Q$ n , 

n 5 = ^ 5 (Q 5 (l) 2 + Q 5 (2) 2 ) +I /t3 _[Q-(q 1 )f 

+ E^5(i)Q 5 (2)(Q3(qi)r> (9) 

n 

where the k's are stiffnesses and we put angled brackets around Qs(qi) to indicate that 
these higher than quadratic terms give rise to effective quadratic terms written here when 
Qs((li) (whose value minimizes the free energy) is regarded as the constant ((^(qi)). To 
make the above Hamiltonian invariant under all the operations of the parent tetragonal 
system, we would have to include analogous terms involving 7^4<5^(qi) = Q% (q 2 ), where 
q 2 — (1/2? — 1/2,0). However, since we are considering a scenario in which Q% (q 2 ) is zero, 
we do not need to write these terms here. Since &i = due to wave vector conservation, Eq. 
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IS 



n 5 = ^Us-hiQzMY' 



[Q B {1) - g 6 (2)] : 



4(K 5 + MQ 3 -( qi )) 2 ) [g5(1)+ 2 Q5(2)]2 
= ^ y r 2 + i/^z 2 , (io) 

where K = [Q 5 (l) - Q 5 (2)]/V2, and Z = [Q 5 (1) + Q B (2)]/V2. This notation (of Table HV]) is 
motivated by the labelling of the axes in the orthorhombic phase (see Eq. (jSJ) and Fig. [3]). We 
therefore see that the presence of a nonzero value of (qi)) lifts the degeneracy in energy 
between the two Q§ modes and leads to the proper linear combinations Y and Z as the new 
normal modes of the Cmcm structure. The analogous splitting occurs in the same way for 
the Q§ modes and leads to the redefinition of normal modes as Y + = [Qe~(l) — Q^(2)]/\/2 
and Z+= [Q+(l) + Q+(2)]/ v / 2. 



1. Induced Absorption Modes 

We next investigate the mixing of modes due to the nonzero value of (Q^(qi))- The 
Hamiltonian we consider for the mixing of absorptive modes is ^3 + + V5, where 

z k 



+ E(Q 3 (qi)) n f« (n) ^ + b (n) v n Y 



+c (n) iy n g 3 -(o)j , (11) 

where U n , V n , W n are single phonon operators (to be determined) which, as we will see, 
become optically active by virtue of their coupling to the optically active Z, Y, or Q^{0) 
modes by the interaction V5 at order (Q3 (qi)) n - To identify these newly active modes, we 
analyze the symmetry of the coupling terms in Eq. ( TTTT) . Note that wave vector conservation 
indicates that terms in Eq. (II ip with n even involve newly active modes at zero wave 
vector and those with n odd involve newly active zero wave vector modes which, before the 
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distortion, were at wave vector q x in the tetragonal Brillouin zone. The character table for 
the X wave vectors was given in Table [Til 

To see more specifically what the Hamiltonian of Eqs. ffTO]) and (fTTj) implies, consider 
only the terms involving U n and Z: 

z z n 

+ E« (n) (^3(qi)) n ^. (12) 

n 

In carrying out perturbation theory relative to the tetragonal phase we use the correct linear 
combinations Z and Y. At the moment we ignore the perturbation in the energy due to this 
interaction which was addressed by Eq. f llOp . Our interest here lies in the mixing of modes. 
When the quadratic Hamiltonian of Eq. ( fl2l) is diagonalized, the new modes U n perturbed 
through their coupling to Z will now be 

(13) 

K>Un ~ K Z 

By analogy we will also have new modes 



c (n) (g 3 (qi)) n ^(o) 



What this means is that the modes U n , V n , and W n (which we will identify in a moment) 
are optically active because they are linear combinations with the optically active modes Z^, 
Yfc, and Qs(0)k, respectively, whose polarizations they inherit. For quantitative work (see 
Appendix A), it is necessary to keep track of the coupling to the different optically active 
modes, Z^, Y^, and Qj(0)fe. For simplicity, since we are mainly interested in which modes 
are newly induced, we need not keep the sum over k. 

Near T>, (Q3(qi)) will be proportional to (T> —T) 13 , where (3 = 1/2 in mean field theory. 
Thus the temperature dependence of terms with different powers of this variable will be 
distinguishably different. Analogous terms for the scenario when the wave vector q 2 has 
condensed can be obtained by applying the operator TZ^ to H3 + H5 + V5. Now we study 
how symmetry restricts the U's, V's and W's of Eq. (fiT|) . Since 7I4 is not a relevant 
symmetry operation, we require ^3 + "Hs + V5 to be invariant under X, m^, and m z . If n is 
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odd (even), then U n , V n , and W n are at wave vector qx (wave vector zero). To be invariant 
under inversion, U n , V n , and W n have to have parity (— l) n+1 under inversion. Also U n and 
V n must be even under m z = (xyz) and W n must be odd under (xyz). Finally we have to 
consider the effect of m d = (yxz). From Tables HI and HT1 we find that 

m d Y = -Y , m d Z = Z , 
mdQz (0) = Qs(0) , m d Qz (q x ) = -Q 3 (q x ) . 

For V5 to be invariant under m^, we therefore must have 

m d U n = (-) n U n , m d V n = (-) n+l V n , 

m d W n = (-) n W n . (14) 

Note that all the above constraints depend on n only through (—1)™, so that to find the 
leading temperature dependences just below T>, we only need to consider two distinct cases, 
n = 1 and n = 2. First consider n = 1. Then 

XU^qx) = C/i(qi) , m^i(qi) = C/i(qi) , 

mdC/^qO = -f/i(qi) , XT4( qi ) = K( qi ) , 

^Vl(qi) = Vi(qi) , m d Vi(qi) = Vi(qi) , 

JWi( qi ) = ^i(qi) , m z W x (ci x ) = -WMi) , 

maW^qi) = -Wi(qi) 



1 



so that 

where a±, /3\, and 71 are constants whose values are not fixed by symmetry. The meaning of 
U\ is that because of the allowed coupling between Z and any operator having the symmetry 
of Q% (qi) this coupling will render all such modes optically active with polarization Z 
inherited from coupling to Z. This coupling can be analoyzed for each Wykoff orbit, as we 
do in Table IVnTl To summarize: the coupling at first order in (Qj(qi)), which we denote 
[and later the superscript 2 indicates second order in (Q3 (qi)}] is 

rn(X+)] 



VP = (Q 3 (qi)> 



k=l 
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n(X+)] 

+ E bjpYQttqj)* 

k=l 

+ E 4 1) Q 3 -(0)Q 4 + (q 1 ) fe 
fc=l 



(15) 



The fact that there are more than one occurrence of (^(qi) affects the intensity but not 
the number or symmetry of new modes, Q 2 (q,i)k, Qi~(qi)fc, and QX (qi)fc- A more complete 
treatment of this coupling is given in Appendix A. 

We now obtain the n[r(qi)] required in Eq. f JT5|) . Using standard methods we find for 
wave vector q 4 that the reducible representations of the displacement vectors over each orbit 
have the decomposition 

r[a(q x )] = X 2 -(q 1 )+X 3 -(q 1 )+X 4 -(q 1 ) , 

r[e( qi )] = X 2 -(q 1 )+X 3 -(q 1 )+X 4 -(q 1 ) 

+X 1 + (q 1 )+X+(q 1 )+X 4 + (q 1 ) 
r[^( qi )] = 2Xr(qi) + 2X 2 -(q 1 )+X 3 -(q 1 ) 

+X 4 -(q 1 ) + 2X 1 + (q 1 ) + 2X 2 + (q 1 ) 

+X 3 + (q 1 )+X 4 + (q 1 ) , (16) 

Therefore, for example, nfX^qi)] assumes the values 0, 1, 1 for orbits a, e, and g, respec- 
tively, leading to the results in the second column of Table IVIII1 

We also show figures of the induced absorption modes in Fig. HI The easiest way to 
generate a mode of a given symmetry (X£ or X^) is to consider each Wyckoff orbit in turn. 
For a given orbit assign one site an arbitrary displacement vector (u x ,u y ,u z ), then generate 
the displacements of the other sites using the characters of Table [III In doing this one will 
find that one or more of the initial displacements can not be nonzero when one includes the 
restriction that the wave function have the given wave vector, in this case qi. 

Now we turn to the case n = 2, for which 

XU 2 (0) = -U 2 (0) , m z U 2 (0) = U 2 (0) , 

m d U 2 (0) = U 2 (0) , 1V 2 (0) = -y 2 (o) , 

m z V 2 (0) = V 2 (0) , m d V 2 {0) = -V 2 (0) , 

XW 2 {0) = -W 2 (0) , m z W 2 (0) = -W 2 (0) , 

m d W 2 (0) = W 2 (0) . 



TABLE VIII: Number of additional optically active modes Q in the orthorhombic Cmcm phase 
for wave vector qi. These are given for the Wyckoff orbits a, e, and g, using the decompositions 
of Eq. © for n\T^ (0)] and n[r+(0)] and of Eq. (USD for n[X+( qi )] and nLY~( qi )]. Pol denotes 
polarization in the notation of Table IVl The order (in the admixture of the wave function) is given 
as proportional to cr n , where a = (Q3 (qi)). If the wave function has an admixture of order a n , 
then the cross section has a contribution of order o 2n . TOT refers to the total for 2 a orbits, 3 e 
orbits and one g orbit. Note that the last row is consistent with the last two rows of Table [VTT1 





Absorption (Odd parity) 


Raman (Even parity) 


Mode 


«#(qi) 


<#(qi) 


Qt(<ll) 


Q 2 (0) 


QF(qi) 


Q 2 (qi) 


Q 3 (qi) 


Qa (qi) 


Q+(0) 


Pol. 


X 


z 


Y 


X 


xz 


XY 


X 2 , Y 2 , Z 2 


YZ 


y2_ Z 2 


Order 


a 


a 


a 


a 2 


a 


a 


a 


a 


a 2 


RP(a) 

















1 


1 


1 





RP(e) 


1 





1 








1 


1 


1 





RP(g) 


1 


2 


2 





2 


2 


1 


1 





TOT 


4 


2 


5 





2 


7 


6 


6 






The solutions for Z7 2 (0) and ^(0) do not produce new modes but instead add to the intensity 
of already existing modes. New modes arise from 1^(0) = Q% (0), so that the relevant second 
order coupling perturbation is 

n[r 2 -(o)] 

K (2) = (Qs (qi)) 2 E 4 2) ^(o)Q 2 -(o) fe , 
fc=i 

where rapT 1 ^ (0)] is given in Eq. (jJJ). We summarize our results in Table rvTTTl Note that all 
the newly induced modes have absorption intensity proportional to (Q$ (qi)) 2 only because 
the low symmetry sites of I4/mmm (which allow the OP Q^(0) to appear) are not occupied 
for CMO or CTO: n[r^(0)] = for CMO or CTO. But an argument equivalent to what we 
have given here would be needed to establish this result. 
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FIG. 4: (Color online) Modes of irrep left and X^~, center, and , right, all at tetragonal 
wave vector qi. These modes become absorption-active in the Cmcm phase. The arrows represent 
displacements in the x-y plane and + and — are the signs of displacements along the z-axis. Here 
the thick-lined squares delineate the equatorial oxygens (at z = ±r) and the uncircled arrows and 
±'s (if any) in the center of the square indicate the displacement of the top apical e-site oxygens 
in the upper layer of the bilayer at z = x an d the circled arrows or ±'s (if any) refer to the bottom 
apical oxygens in the lower layer of the bilayer at z = All the e-site ions move along the same 
axis as the apical oxygens, but with independent amplitudes. There are two modes of symmetry 
corresponding to independently choosing the amplitudes of the large black and small red arrows. 
If the large black arrow is dominant, this mode is the (0, 0, 0) rotation about the tetragonal axis. 
There are five modes because the displacements of the apical oxygen, the other two e-site ions 
(not shown), and the two components of the equatorial oxygens are all independent parameters. 
There are four Xf modes because the tilting, the displacements of the apical oxygens and of the 
two other e-site ions are independent parameters. (Some of the modes, such as Q~± (qi) and (q.4), 
which grossly distort the octahedra will have relatively high frequencies.) Without the admixtures 
of order (Q3 (qi)) all these modes are of the wrong symmetry to be optically active. However, the 
admixing, e. g. Eq. (|15p in the Cmcm phase makes these modes optically active. 
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2. Induced Raman Modes 



Here we consider the admixing of the Raman active modes (see Tabled]) from irreps r+, 



with n = 1,2,4,5: 



Vc = E(^3 (qi) ? 



>Z + T n + b n Qt(0)U n 



+c n Q+(0)V n + d n Qf(0)W n + e n Y + S n 

where T n , U n , V n , W n , and S n are the modes which are Raman active by virtue of their 
coupling proportional to (Qs)( c li) n to modes of the tetragonal system which are Raman 
active. For n = 1 invariance with respect to X, m^, and m z leads to 



So 





= -5i(qi) , 


^d5"i(qi) = S 1 (q 1 ) , 


rrizS^qi) 


= -^i(qi) , 


XT 1 (q 1 ) = -T 1 (q 1 ) , 




= -^(qx) , 


m 2 Ti(qi) = -Ti(qi) 


^(qx) 


= -Ux(qx) , 


m rf C/i(qi) = -C/i(qi) 


m z Ui(qx) 


= tMoa) , 


XT4( qi ) = -^(qj) , 




= ^( ql ) , 


m z Vi{(\i) = Vi(qi) , 


2Wi((u) 


= -^i(qi) 


JTidWi(qi) = -Wi(qi 


^Wi(qi) 


= Wi(qi) • 





^i(qi) = Qi (qi) , ^i(qi) = Q x (qi) , t/i(qi) = Q 3 (qi) 
K(qi) = Qr(qi) , Wi(qi) = Qi-(qi) . 



Thus 



(n=l) 



C 



(Oi"(qi)> 

n(X 3 ") 



n(X-) 



E ^^rc*: 



fc=i 



+ E ^Q 4 + (o)Q 3 -(qi)* 
jt=i 

"(*4~) 
fe=l 



+ E c^WQKqa)* 



21 



n(X-) 

k=l 

nX~) 

+ e«y+g 2 -( qi ) fe 
fc=i 



where the sums over the number of appearances of the Raman active tetragonal modes is 
omitted, in analogy with Eq. ([15]) . The fact that (0) and Qt(0) are both coupled to 
X^(qi) k means that each X 3k mode has two independent admixtures, one, ^4(0), giving a 
contribution to the Raman tensor a(ij + ji) = a(KK — J J) and the other, Qf(0), giving a 
contribution to the Raman tensor bkk + c(ii + jj) = bit + c( J J + KK). As a result, the 
X^ k {c\x) modes have a Raman tensor bll + (c — a)(JJ) + (c + a)KK, so that X^ k (qi) has 
a Raman polarization we denote (X 2 , Y 2 , Z 2 ). 

For n = 2 invariance with respect to X, m^, and m 2 leads to 

XS 2 (0) = 5 2 (0) , m d S 2 (0) = -S 2 (0) , 

m a 5- 2 (0) = -S 2 (0) , XT 2 (0) = T 2 (0) , 

m d T 2 (0) = T 2 (0) , m,T 2 (0) = -T 2 (0) 

XU 2 (0) = U 2 (0) , m d U 2 (0) = U 2 (0) , 

m z U 2 (0) = U 2 (0) , 1V 2 (0) = V 2 (0) , 

m d V 2 (0) = -V 2 (0) , m z \/ 2 (0) = V 2 (0) , 

XW 2 (0) = W 2 (0) , m d W 2 (0) = Vr 2 (0) , 



m z W 2 (0) = W 2 {0) . 



So 



s 2 = 


Y+(0), T 2 = Z + , 


u 2 = 


aQt(0)+PQt(0) , 


v 2 = 


a'Qt(0)+P'Qt(0) 


w 2 = 


a"Qt(0) + /3"Qi(0) 



The only new irrep here is T^O) with OP Q3 (0). So we have the results of Table IVIIII 
As was the case for the absorption cross section, all the newly induced modes have their 
Raman scattering cross section proportional to (Qj (qi)) 2 only because the low symmetry 
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sites of I4/mmm (which allow the OP Q£(0) to appear) are not occupied for CMO or CTO. 
We emphasize that this result is not a general result for the space group I4/mmm distorting 
into space group Cmcm. 

In Figs. [5] and |6] we show the new Raman modes in the Cmcm phase. 




£ x = -l 1 2 x = -10 1 2 



X l X 2 

FIG. 5: (Color online) As Fig. |4] except for modes of irrep X-f, left and X^T, right, both at 
tetragonal wave vector qi . These are the modes which become Raman active (see Table IVIIIj) 
due to admixtures of tetragonal Raman active zone center modes. In mode the sites on the 
z axis do not move, so there are two modes of this symmetry corresponding to the choice of the 
amplitudes of the x and y displacements. There are 7 modes of symmetry X^(qi) because there 
are two ways to choose the amplitude of the arrows (the displacements of the equatorial oxygens) 
and one way to independently choose the amplitudes of the displacements along the z-axis for each 
of the 2 a orbits and for each of the 3 e orbits of such sites. These breathing modes grossly distort 
the octahedra and therefore have relatively high frequency. 



3. Temperature Dependence of Mode Energies and Absorption Intensities 

Near the phase transition at temperature T = T> where the order parameter (Qsi^i)) 
becomes nonzero one has 

u n (T) = J°\T) + 6ou n (T) , 

In(T) = I$»{T) + 6I n {T), (17) 
23 



FIG. 6: (Color online) As Fig. [5] except for modes of irrep X% , left and X± , right, both at 
tetragonal wave vector qi. The equatorial oxygens (g orbits) tilt around the [110] direction for 
irrep X^(qi) and around the [110] direction for X^(qi), in each case with a single amplitude. 
The other ions on the z-axis move along the [110] axis for X^(qx) and along the [110] axis for 
X^iqi), in each case with an independent amplitude. Thus for both irreps we have 6 modes (one 



for each orbit). In the notation of Refs. 



20 Q 3 (qi) is a (<&, $,0) rotation and Q 4 (qi) 
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is a ($,—$,0) rotation. The wave vector regulates how one bilayer is structured relative to the 
adjacent bilayer. So these two modes have slightly different frequencies because although their 
intrabilayer interaction energies are the same, their interbilayer interaction energies are different. 



where the first term is the smooth temperature dependence and the second term is the 
anomalous term at the phase transition which is only nonzero for T < T> and will involve 
powers of (Q^ (qi)) = & ~ |2~> — T\ @, where the order parameter exponent /3 is 1/2 in mean 
field theory, but fluctuations will cause it to be less than 1/2. 

First consider the mode energies. For all the nondegenerate q = modes, the perturbation 
of the form given in Eq. (ITTT) . will lead to an energy shift from terms with n — 1, which, in 
second order perturbation theory, gives a contribution to 5u(T) of order a 2 ~ \T > — T\ 2 ^ . 
For the degenerate r^(0) modes, 5ui(T) has a contribution from first order perturbation 
theory, but in that case, since b\ of Eq. (jHJ) vanishes, the contribution to 5u(T) is again of 
order a 2 ~ |T> — T\ 213 . These results are shown in the left-hand panels of Fig. [71 It should 
also be noted that the mode energies can be modified by quartic terms in the Hamiltonian 
which we have not so far considered and which are of the form 
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35] 



V 4 = Y. d nUlQ^x) 
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which will also lead to a contribution to 5uj(T) of order a 2 ~ |T> — T\ 2/3 . It is plausible, 
but not a certainty, that the sign of the anomalous contribution would cause a stiffening 
of the modes as the temperature is reduced as is depicted in lower left-hand panel of Fig. 
CD If V4 were to be neglected, then one would have a direct correlation between SI n (T) and 

<MT).|afi| 

Now consider the mode intensities in the presence of admixtures written in Eq. (1131) . 
The standard scenario in perturbation theory is that a perturbation in first order introduces 
new components into the wave function without changing the unperturbed component of 
the wave function. In second order perturbation theory the amplitude of the unperturbed 
component of the wave function can be modified. Here, this tells us that if the mode U n is 
already dipole allowed, then the intensity will be modified at order a 2 . If the mode U n is 
not allowed in the absence of the perturbation then its intensity will be proportional to the 
square of the amplitude of the admixture. Furthermore, new modes that appear at first order 
ina will have intensity of the form of Eq. ([IT]) with/^(T) = and 6I(T) ~ a 2 ~ |T>-T| 2/3 . 
New modes that only appear at second order in a will have intensity of the form of Eq. f fT7|) 
with I^\T) = and 51 (T) ~ cx 4 ~ |T> - T| 4 ^. From Table EM we see that for CMO 
or CTO there are no such modes that only appear at order a 2 . This is a special property 
of the RP system. Had we had a more general tetragonal system with low symmetry 
sites, then there would be modes which only appear at order a 2 and which would have 
51 (T) ~ er 4 ~ \T > — T| 4/3 . Our results are summarized in the right-hand panels of Fig. [7J 



IV. CMC2i PHASE 



A. Analysis based on Cmc2i symmetry 

Since Cmc2i does not have a center of inversion symmetry, parity is not a good quantum 
number and Raman and absorption modes are mixed. In Cmc2i the general orbit is half as 
large as for Cmcm. However, the unit cell then contains twice as many orbits. The sites in 
Cmc2! are given in Table IIXI and the character table for q = is given in Table |2D 

Using the data in Table |X] we find that 

r(a Cmc2l ) = 2r x + 2r 2 + r 3 + r 4 , 
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FIG. 7: (Color online) RP modes. Left panels: mode energy versus temperature, where T> is the 
temperature at which (Q3 (qi)) becomes nonzero. Right panels: the mode intensity versus temper- 
ature. Upper left: removal of degeneracy of the modes. Lower left: temperature dependence 
of a nondegenerate mode energy. Upper right: temperature dependence of the intensity of a non- 
degenerate mode which is absorption allowed in the tetragonal phase. Lower right: temperature 
dependence of nondegenerate mode intensity for modes which have zero intensity in the tetragonal 
phase. In all cases, the discontinuity in slope is drawn assuming the mean-field result, (3 = 1/2, 
for the order parameter exponent. Otherwise, the lower temperature additional contribution is of 
order |T> — T| 2 ^. 



TABLE IX: Operators for Cmc2i. 

E = (x,y,z) 2 z = (x,y,l/2 + z) 
m x = (x,y,z) m y = (x,y, 1/2 + z) 

4 4 

r(e CmC 2i) = 3 r n , r(g Cmc2l ) = 6 ^ r n , (18) 

n=l n=l 

and the dipole moments and Raman tensors of each irrep are given in Table IS Note that 
once we have identified the absorption modes, the identification of the Raman active modes 
follows, except for those of irrep IV Using Eq. (fT8|) we find the number of optically active 
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TABLE X: Character table for the little group at q = for the Cmc2i lattice. The last three 
columns are the characters for the various tetragonal orbits on the Cmc2i structure. Orthorhombic 
coordinates are used throughout this table. P gives the nonzero component of the dipole moment 
operators for the irrep in question. The Raman tensors R are R(Ti) = dll + eJJ + fKK, 
R(r 2 ) = g(JK + KJ), R(r 3 ) = h{IJ + JI), and R(r 4 ) = i(IK + KI). 



Operator 


Ti 


r 2 


r 3 


r 4 


a 


e 


g 


(X,Y,Z) 


1 


i 


i 


i 


6 


12 


24 


(X,Y,Z) 


1 


i 


-l 


-l 


2 








(X,Y,l/2 + Z) 


1 


-l 


i 


-l 











(X,Y,l/2 + Z) 


1 


-l 


-l 


i 











P 


aK 


bJ 


None 


ci 









modes in the Cmc2 1 phase given in Table IXII 

B. Analysis based on lowering of Cmcm symmetry 

Now we perform the analysis of mode symmetry using the mixing of modes in terms of 
the OP's present in the Cmc2! phase, namely Q% (qi), Z = [T^ + T^ 2 ]/V^, and Qt(<ii), 
where we assume the OP's of the family of qi. (Afterwards the analogous terms associated 
with q 2 can be obtained by applying the four-fold rotation TZ^ to the coupling we construct 
for qi.) At the lower transition we can mix using all the condensed order parameters, 
Q3 (qi)) an d Z. The first question is: can we induce Y = <3^i(0) — <3^2(0)? The 

answer is no. Note that cr^Z = Z and cr^Y = —Y. Also (TdQ^i^li) — ~ (T dQ^{ ( Ai) and 
(TdQt (Qi) = ~Qt( c li)- But for wave vector conservation we must have a total even number 
of powers of these two operators at qi. Hence there can be no term in the free energy which 
is linear in Y. So Y is not induced by the other order parameters. 
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TABLE XI: As Table IXl Optically active modes of the Cmc2i phase. TOT gives the total number 
of modes of each symmetry found using Eq. (|18p and Table [Xj We copy the number of modes in 
the Cmcm phase from Table IVTI1 and Addl gives the number of new modes which appear when the 
Cmc2i phase is reached, as given in Eqs. (|20p and (|24p . Here R refers to Raman active modes and 
P to dipole-active modes. This entry includes only the additional modes we explicitly calculated 
for the dipole active modes and for the Raman active modes which are not dipole active. Since the 
total number of modes is given for both Cmcm and Cmc2i, one can deduce that, for instance, the 
number of additional Raman active modes with Raman tensor i{IK + KI) is 17-7=10. 
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r 3 


r 4 
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aK 


bJ 


None 


d 


TOT 


19 


19 


17 


17 


Cmcm 


10(R), 9(P) 


7(R), 12(P) 


12(R) 


7(R), 10(P) 


Addl 


10(P) 


7(P) 


5(R) 


7(P) 



1. Induced Absorption Modes 

So optically active absorption modes are admixed by the perturbation 

v{m)= E(^(qi)) fe (^(qi))'(^) m 

k,l,m 

x {^klmUklmZ + b k i m V k i m Y + CklmWklmQz (0)1 , 

where again the U's, Vs, and W's are the single phonon operators which will become 
optically active due to this mechanism and which we wish to identify. We insert the angle 
brackets to emphasize that these higher order terms give rise to an effective bilinear coupling 
of the type that mixes modes. In principle any perturbation has to be invariant under all 
the symmetries of the tetragonal phase. However, applying a four-fold rotation takes qi into 
q2 (whose order parameters are zero) and {Z) into (Y) (which is also zero). Thus we require 
that V(qi) be invariant under only m^, m z , and I, in which case 

TndUkim = (—) k+l+1 Ukim , rridVkim = {—) k+l Vkim , 

"mdWklm = ( — ) k+l+1 Wklm , rn z Uklm = Uklm , 
^zVklm = Vklm , ^z^klm = — ^klm 
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ZUklm — (—) +m+1 Ukl m , IVklm — (~ ) +m+1 Vkl m , 

1W klm = (-) k+m+1 W klm . 



TABLE XII: Independent sets of indices for the operators U, V, or W (which depend on (— ) fc + m 
and (— ) k+l ). "Cmcm" indicates equality to cases previously considered for Cmcm. The equal sign 
indicates equality in view of (Z) oc (Q% (qi))(Q^(qi))- Cases which lead to identical operators 
with different prefactors are indicated by "oc." Cases which produce "new" operators are indicated 
by *. 
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We use the result of Ref. |9j that, due to the cubic coupling between the order parameters, 
one has 



(Z) = a(Q+( qi ))(Q 3 "(qi)) 



(19) 



where a is a constant. Note the following: a) increasing k, I, or m by two does not induce an 
additional operator and b) solutions for U, V, and W, depend on k + l and k + m. Therefore 
we have the results of Table IXLTI of which we only consider cases A, B, and E. We find that 

tfooi = cQ+(0) + dQt(0) , V 001 = c'Q+(0) + d'Qt(0) , 
Wooi = Z + , Uxoi = C/qio = QaM , 

V im = V 010 = Qj(qi) , W 101 = W Q10 = Qr(qi) • (20) 
From these three cases we get (for wave vector qi) 

V(qx) = (Q 2 + (qi))(Q3(qi)){«iQ 3 (0)^ + + z[ 



a 2 Qt(tt) + a 3 Qt(0)} + Y[a 4 Qj{0) + a 5 Qj 
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+(QtM)\a*ZQ3 (qi) + arYQlM 

+a 8 Q 3 -(o)Qr(qi)} + (QtMKQs (qi)> 2 { 

a 9 ZQ 3 _ (q 1 ) + ai FQ 4 (qi) + flnQg (O)Qi (qi) J 

By applying 7^ 4 to this we get 

V(q 2 ) = (QtM}(Q 3 (cL2)}^a 1 Q 3 (0)Y + + Y[ 

a 2 Qt(0) - a 3 Q+(0)] + Z[a 4 Q 2 + (0) - a 5 Q+(0)] 

+(Qt(<i2)){aeYQz (q 2 ) + a 7 ZQj (q 2 ) 

+a 8 g 3 -(o)Qr(q 2 )} + (Q^(qa))(Qj (q 2 )) 2 | 

agFQg (q 2 ) + ai ZQ4 (q 2 ) + anQ^ (0)Q^ (q 2 ) \ 



(21) 



(22) 



<?. Induced Raman Modes 

Finally, we carry out the same analysis to determine the new Raman modes of polarization 
X Y , since these modes are not dipole active and hence are not accessible to the previous 
calculation. To get polarization X Y = z t (x t — y t ) we must induce the mode Y + of irrep 
I5 . We use the results of Table Kill to restrict our attention to the coupling for the cases 
k = I = m — 1 = (together with k = l = m+l = l) and k = I — 1 = m = (together with 
k — l + l — m — 1) so that 

V = R.Y^aiZ) + /3(Q^(qi))(Q 3 _ (qi))] 

+R 2 Y + [ 1 (Qt( qi )} +8(Qz(c ll ))(Z)] , (23) 

where Rk are the single phonon operators which we determine by requiring that V be an 
invariant under I, rrid, and m z . we find that 

R x = Qi(0) or £4(0) , R 2 = Q+(q x ) . (24) 
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TABLE XIII: Additional optically active modes (which we refer to as 'new') which appear in the 
Cmc2i phase for wave vector qi from Eqs. (|20p and (|24p . The values of n[(r(0)] for each orbit are 
given in Eq. ([3]) and those of ra[(T(q)] are given in Eq. f)16[) . TOT is the total number of modes 
from 2a orbits, 3e orbits, and one g orbit. The values of P (and the associated values of the Raman 
tensor R) are from Z, Y and Q^(0) in Eq. (|2ip . The values of the Raman tensor R for Q4 (0) 
is inherited from Y + in Eq. (|23p . In the last row the entry 1 indicates that the Raman tensor is 
K = ail + bJJ + cKK. 
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1 


1 
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YZ 


XZ none XY 


1 


YZ 


XZ 
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The number of modes we admix into R\ is equal to the number of times either IY(0) or 
rj(0) appears, which from Eq. (J3]) is one for each g orbit or once in all. The number of 
modes we admix into R2 is equal to the number of times irrep ^(qi) appears, which from 
Eq. (j!6p is one for each e orbit and one for each g orbit, or four in all. So there are five 
additional X0Y0 Raman modes induced in Cmc2i. 

We could give figures showing the displacements in the new optically active modes which 
appear in the Cmc2i phase. Actually, the new modes at wave vector qi have already been 
shown in Figs. HJ El and El but now many of these modes are both absorption active and 
Raman active. We leave it to the reader to draw the new modes at zero wave vector using 
the same technique as used above. 
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3. Summary of Optically Active Phonons. 



In Table IXIVI we list the irreps of the tetragonal phase which are active or become active 
due to the condensation of order parameters (q), Qt(<l), and r^~(0). However, one can 
show that the intensities for wave vector q 2 are obtained from those for wave vector q x by 
simply replacing qi by q 2 . 

TABLE XIV: The irrep labels are according to the tetragonal structure. "Pol" denotes polariza- 
tion, where capitals (lower case) refers to orthorhombic (tetragonal) coordinates, "Int" denotes 
intensity, and "Modes" denotes the number of modes. The columns labeled q n apply when the 
wave vector q n is selected, oi = E/ /3fc;(Q^(q)) 2 , <?2 = E/, m 7fc/m(Q^(q)) 2 (Q^ m (q)) 2 ' and 
a 3,k = Y.im(Qt,i( c i))' 2 lkim + $kim {Q31 (q)) 2 ] 2 , where k is the mode index and a, /3, and 7 do not 
depend on the OP's. 
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V. MAGNETOELASTIC MODES 



Here we investigate the mixing of dipole active phonons with magnons for the RP 
system containing the magnetic ion Mn 4+ . The motivation for this investigation is that 
photon absorption by magnons normally proceeds via a magnetic dipole matrix element, 
whereas absorption by phonons proceeds via the very much larger electric dipole matrix ele- 
ment. In magnetic ferroelectrics, magnons into which phonons are mixed have been dubbed 
"electromagnons." 20[ Such magnons not only have a much enhanced absorption cross sec- 
tion, but, because their energy is often uruch lower then the phouou energies, they can also 
lead to large anomalies in the static dielectric constant. [23] Potentially the same effect is 
possible here if the magnons have a nontrivial coupling to phonons. 

To start we review what is known about the ma gne tic structure of CMO which appears 
below the magnetic ordering temperature of 115K.[37| The OP's are[37| Gx(qi), CV(qi), 
and Fz, or equivalently {9] Gx{<l2), Fy, and Cz(c^), where G, F and C are the Wollan- 



Koehler magnetic structural descriptors 



38| for the magnetism of a single bilayer: F denotes 



a ferromagnetic structure, G and C are the antiferromagnetic structures illustrated in Fig. 
IU The wave vector distinguishes between the two possible stacking (of G and C) as one 
goes from one bilayer to the next,|9 ; ] as shown in Fig. El As mentioned below Eq. (CQ) one 
can select the wave vector by cooling through T> in the presence of a shear stress. 




FIG. 8: Spin states of a bilayer of Mn ions at z = ±£. The plus and minus signs represent the 
signs of any component of spin. Left: the "G" configuration. Right the "C" configuration. The 
uncircled symbols are in the planes at z = ±£ and the circled symbols represent the spins in planes 
at z = 1/2 ±£ when q = qi- For q = q2 the circled + and — signs are interchanged. The tetragonal 
axes are shown at right. 



The plan of our analysis is as follows. We will construct the most general trilinear 
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coupling which involves two magnetic OP's (in order to satisfy time- reversal invariance) and 
one phonon operator whose wave vector is either at the zone center or at an X point. We 
will require invariance with respect to all the symmetries of the tetragonal lattice. (As a 
short-cut we fix the wave vector to be qi. The terms involving q 2 can be determined later by 
applying a four-fold rotation to those we find for qi.) We then develop an effective bilinear 
coupling by replacing one of the magnetic OP's by its value at the minimum of the free 
energy. 

Thus we analyze in detail the cubic interaction (involving the wave vector qi) 

VME(qi) = Y^[^ c <*M m pQ*M 

+B aP7 C a (q 1 )Np(q 1 )Q-(0) 
+C a ^M a Np{q 1 )Q-{q 1 )} 

k=l 



(25) 



where M denotes F, N(q) denotes G(q), and we chose the arguments of the phonon op- 
erators so as to conserve wave vector. We also used the transformation properties of the 
magnetic OP's under inversion (see Table IXVl) to infer the inversion signature of the Q's in 
Eq. ([25]). 



TABLE XV: The symmetry for components X, Y, and Z of the magnetic OP's (which are pseu- 
dovectors). Here Sx = S z , Sy = S x — S y , Sz = S x + S y , where capitals refer to orthorhombic and 
lower case to tetragonal coordinates. The results do not depend on the choice of wave vector. In 
the text we use the standard notation in which the net magnetization, here F, is denoted M and 
the staggered magnetization, here G(q), is denoted N(q). 
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To get the effective coupling between magnons and phonons we only analyze special terms 
of Eq. (|25p . We take the expectation value (that minimizes the free energy) of one of the 
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TABLE XVI: The symmetry of the magnetoeleastic trilinear interaction for wave vector qi. All 
components are with respect to orthorhombic coordinates. The variables C and N are at wave 
vector qi, as are the irreps X/~. Under (m z ) we give the effect of (m z ) on each of the two 
magnetic OP's. Then the irrep(s) is identified so as to make the trilinear interaction invariant 
under and m z and in the next-to-last column we identify the polarization P of the absorption 
mode (if any). Only irreps that are optically active are included in the analysis. In the last column 
we give the order of magnitude of the term, assuming the staggered moment to be unity. 
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magnetic OP's and therefore arrive at a quadratic interaction between the other (active) 
magnetic OP and a phonon OP. To get a magnon-phonon coupling out of this we select the 
component of the active magnetic OP to be transverse to its equilibrium value. Since the 
nonzero equilibrium component of N(qi) is its X component, the transverse components 
of N(qi) are its Y and Z components. These are the components of N(qi) which, in the 
Holstein-Primakoff formulation, create spin waves. The operator Q creates or destro ys a 
phonon. These terms which mix magnons and phonons give rise to "electromagnons." 20- 
231 ] We treat the other trilinear terms similarly, keeping only those in which one magnetic 
OP is longitudinal and the other magnetic OP is transverse. Thereby we have the cases 
listed in Table IX VII for wave vector q x . The results for wave vector q 2 are quite similar: 
the polarizations Y and Z are interchanged and in the intensities qi is replaced by q 2 , as in 
Table IXIVI It is relevant to try to estimate which of these terms are the most significant. 
For that purpose we keep in mind the fact that the OP's M and C(qi) are roughly one tenth 



35 



as large as N(qi). 37]] In addition the optical activity of the phonon depends on whether it is 
allowed in the tetragonal phase or depends, for its activity, on the distortions of irreps X% or 
X% (in this connection we ignore Cf). It seems likely, then, that the largest magnon-phonon 
coupling will involve coupling the transverse moment of C(qi) to either X or Z polarized 
zone center phonons. But for completeness, we list the other possibilities. Table IXVII only 
gives results for the case when wave vector q! is selected. Results for wave vector q 2 can be 
obtained by replacing Z polarization by Y and qi by q 2 , as is done in Table IXlVl 

These couplings can be detected either by anomalies in the static dielectric constant 
when magnetic order appears, 23j or by anomalous electric dipole intensity in magnons due 



to their coupling with optical phonons. 
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VI. CONCLUSIONS 



In this paper the effect on the optical properties of zone center phonons of the sequen- 
tial lowering of symmetry by structural distortions from the high symmetry tetragonal 
(I4/mmm) phase has been studied. In addition, the coupling of the magnons which ap- 
pear at low temperatures in the magnetically order phase of Ca3Mn207 and the phonons is 
explored to identify possible giant enhancements in the magnon cross section in the electro- 
magnetic spectrum. 

Our specific results include 
• 1. We have tabulated the zone center phonon modes with their polarizations which appear 
in the absorption spectrum and in the Raman scattering spectrum of the I4/mmm, Cmcm, 
and Cmc2i phases for 0^3X207 systems, where X is Mn or Ti. 

•2. For modes which become optically active by virtue of a structural distortion, we indicate 
the dependence of their cross sections on the newly emerging order parameters near the phase 
boundary where the structural symmetry is lowered. We also discuss the anomalies in the 
phonon frequencies as the phase boundaries are crossed. 

•3. We discuss the symmetry of the magnon phonon interaction which can lead magnons 
being electric dipole active and hence having anomalously large absorption cross section and 
leading to analogous large magnetic contributions to the dielectric constant. 
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Appendix A: Intensity Calculation 



1. Tetragonal Phase 

In this appendix I discuss in detail how the intensities of the modes are calculated as- 
suming the perturbative description relative to the tetragonal phase is justified. In the 
tetragonal phase we have the optically active phonons of irreps and , governed by the 
Hamiltonian 

i "(17) i n(r~) 

^tet = - E K ± , k (Y k 2 + z 2 k ) + - E «n,*** , 
z fc=i z k=i 

where it is conveneient to work in orthorhombic coordinates, so that modes X, Y, and Z 
are polarized along these respective orthorhombic direction. X k = Qs{0)k, Y k — (Qs — 
Q§ (2)fc)/y / 2, and Z k = (Q^(l) k + Q^(2) k )/y/2 are the order parameters for the modes of 
the optically active and irreps. To be specific: we number the modes in order of 
incresing energy. 

The coupling to photons is via the dipole moment operator which we write as 

n(r-) n(T-) 

P = E P±,k[YkY + Z k Z}+ £ P\\,kX k X , 

fc=l k=l 

where X is a unit vector in the orthorhombic X direction and similarly for Y and Z. So 
the absorption intensity in the tetragonal phase of the kth doubly degenerate Y k -Z k mode is 
Ap\ k and that of the fcth X mode is Ap?, k , where here and below we set A = 1. Since the 
dipole moment matrix elements are subject to a first principles calculation, it makes sense 
to give explicit results in terms of calculable quantities. 



2. Cmcm Phase 



The Hamiltonian that describes the mixing of new optically active phonons in the Cmcm 
phase is 

1 n(X+) 1 n(X+) 

^ Cmcm = - E 4 Xl V(qi)^ E ^QiMl 
z k=i z k=i 

+\ E - { k Xt) Q + Mi)l 
z k=i 
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rn[X+] n[r~] 



n[x+] n[r-; 



L fe=l i=l 



+ E E ftSJ/nof (qi)i 

fc=l i=l 

n[X+] n[r-] 

E E 4^3-(o),Q 4 + (qi)^ 
fc=i i=i 

Q^(qi) is replaced by (Q% (qi)) leading to a bilinear interaction which mixes new modes. 

We will analyze the Z-component of the dipole moment operator. The other components 
are treated analagously. Due to the mixing of modes, the bare mode operator Z k is related 
to the true mode operators (indicated by tildes) 



g* + £ °^ (qi) ' «g(qi)>- 



I K t " - Kx,k 

Thus the Z-component of the dipole moment operator which creates modes Q2 (qi)z denoted 
Pz,q+,i is 

pz,q+,i - p±,k—^+j (<^ 3 (qij> ) ( A1 ) 

and the intensity of this induced mode is \p Z q+ J 2 . The quantities in Eq. (lAip are amenable 
to a first principles frozen-phonon calculation. 
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